Entanglement Entropy and Mutual Information in Bose-Einstein Condensates 
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In this paper we study the entanglement properties of free non-relativistic Bose gases. At zero 
temperature, we calculate the bipartite block entanglement entropy of the system, and find that 
it diverges logarithmically with the particle number in the subsystem. For finite temperatures, 
we study the mutual information between the two blocks. We first analytically study an infinite- 
range hopping model, then numerically study a set of long-range hopping models in one-dimension 
that exhibit Bose-Einstein condensation. In both cases we find that a Bose-Einstein condensate, 
if present, makes a divergent contribution to the mutual information which is proportional to the 
logarithm of the number of particles in the condensate in the subsystem. The prefactor of the 
logarithmic divergent term is model dependent. 

PACS numbers: 03.75.Gs 



INTRODUCTION 



CO . Entanglement, as measured by, e.g., bipartite block entanglement entropy, is playing an increasingly important role 
in the study of condensed-matter or quantum many-body physics, both conceptually and quantitatively. It has been 
used as a very useful and in some cases indispensable way to characterize phases and phase transitions, especially for 
O ^1 phases and quantum phase transitions in strongly correlated fermionic or spin systems (for a review, see Ref. [1]). 

For bosonic systems, studies of entanglement entropy have mostly focused on relativistic free bosonic field theories 
^ . [2-4], which are equivalent to coupled harmonic oscillator systems (for reviews, see Refs. [5] and [6]). 

' In this paper we study the entanglement properties of free non-relativistic Bose gases. In addition to interest in its 
^ own right, our motivation also comes in part from the following consideration. In recent studies it has been shown that 
entanglement is enhanced at quantum critical points [7] and strongly correlated phases with topological order [8] , in the 
form of either violation of area law [2, 3, 9-20], or subleading corrections to the area law that diverges with block size 
CsJ ■ [21-23] (usually in a logarithmic fashion). On the other hand there have been relatively few studies of the behavior of 
I entanglement entropy in states with traditional long-range order [24-26]. In a recent work [27], we calculated the block 
. entanglement entropy of some exactly soluble spin models that exhibit ferromagnetic or antiferromagnetic long-range 
order in the ground state, and found that such conventional orders also lead to logarithmically divergent contribution 
to the entropy. Bose-Einstein condensation (BEG) is perhaps the simplest example of conventional ordering. It is 
thus natural to study its entanglement properties. As we are going to show, a Bose-Einstein condensate (referred to 
as a condensate from now on) indeed makes a logarithmically divergent contribution to the entropy as well. 

Besides the entanglement entropy of the ground state, the entanglement properties of system at finite temperature 
[ are also of great interest. However, the entanglement entropy is only well-defined for a pure state. For a system that 
is described by a mixed density matrix, the von Neumann entropy of the reduced density matrix becomes different 
for the two parts of the bipartite systems. In such cases, there is a natural extension of the entanglement entropy 
that one can work with - the mutual information [28, 29]. We will show that a condensate, when present, makes a 
^ ' logarithmically divergent contribution to the mutual information. 
5^ [ This paper is organized as follows. In Sec. II we study the ground state entanglement entropy of a generic free 
boson model that is translationally invariant [30]. In Sec. Ill we introduce an infinite-range hopping model for bosons 
which is exactly solvable, and calculate the mutual information analytically. In Sec. IV we introduce a long-range 
hopping model for bosons in one-dimension (ID) which exhibits a finite temperature BEG for a certain parameter 
range, then we present a numerical study of the mutual information for this model. In the end, we summarize and 
discuss the results of this paper in Sec. V. 

2. ZERO TEMPERATURE: ENTANGLEMENT ENTROPY OF FREE BOSONS 

Consider a general Hamiltonian of free bosons hopping on a lattice of size L: 

H = - ^Uja\a,, (1) 
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where tij > 0, ai(aj)'s are the bosonic annihilation (creation) operators. If the system is translationally invariant, 
tij = ti^j, then the Hamiltonian can be diagonaUzed by Fourier transformation: 



(2) 



where bk = ^ ^ is the annihilation operator in k space. In most generic cases, the ground state is the k = 

state. At zero temperature, all particles fall into the ground state. For a system containing N particles, the ground 
state is given by: 



(3) 



To consider its bipartite block entanglement entropy, we divide the system of size L in two parts, and label them 
A and B respectively. Let the sizes of each part be La and Lb-, La + Lb = L, and define 



Eat, = at 



Then we can write l^'n) as: 
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where A, = L-^^^L^L^"', \l)^ = ^aJ^lO)^, \N I) ^ = ^^^^^-'jO)^, and |0) = |0)^ ® \0) ^. 
This is an explicit Schmidt decomposition, and therefore the entanglement entropy is readily given by: 



S = -EA(lnAi. (6) 

We are interested in the asymptotic behavior in two limiting cases: (1) the equal partition case; (2) size of B is 
substantially larger than A, i.e.. Lb ^ La- 
(i) Equal partition. La = Lb = ■§■: 



' l\{N-l)\ L^ l\{N~l)\2N + 



Let X ~ — then x G [—"t-, "y], and we can denote A; as A^, = '^,J^.l ,n__ \^riv i h which can be approximated by 

_2x^ 

a Gaussian distribution factor A^; ~ e " when N is large. In the limit N ^ oo, the summation over n (or x) can 
be approximated by an integral. Also in this limit, the Gaussian factor is sharply peaked around .t = 0, the integral 
region can be extended to from minus infinity to infinity. Using the fact that Aa; ~ \{x)dx = 1, we arrive at 



The entanglement entropy is then 



A(x)^W-^e-^. (8) 



1 / Nil \ 1 

Ec,-l \{x)\n\{x)dx = -[l + \n{ — )]=-\nN + 0{l). (9) 
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(ii) Unequal partition. Lb S> La'- 
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If Lb 3> La, i oo, but keep ^ — !■ (n) (fixed), the distribution of A/ approaches a Poisson distribution: 

^' = n(7v-0! ' n ■ ^ ^ 

The entropy of the Poisson distribution, which in this case is our entanglement entropy, is known to be: 

E=l-[l + \n{2TTLA{n))]- —-^ + 0{ ^ ' 



2' ' 12LA{n) \LA{n)Y' ^^^^ 

= 1 [1 + ln(2.A..)] - ^ + O(^) = i IniV. + 

where Na = LA{n) is the average particle number in subsystem A. 

Therefore, we find, in both cases, that the leading term of the mutual information goes as ^ IniV^ for La < Lb- 



3. MUTUAL INFORMATION: ANALYTIC STUDY OF AN INFINITE-RANGE HOPPING MODEL 

In this section, we will study the natural generalization of entanglement entropy at finite temperature: the mutual 
information, which is defined as 

Em = ^{Ea+Eb-S), (12) 

where Ea and Eb are the von Neumann entropy of the reduced density matrices of subsystems A and B, respectively, 
and S is the entropy of the whole system. Note that at finite temperature Ea and Eb are no longer the same due 
to the fact that the system is described by a mixed density matrix. We must emphasize here that our definition of 
mutual information differs from its usual definition [31] by a factor of 2 so that it will converge to the entanglement 
entropy when the system approaches a pure state. 



3.1. Model, spectrum and thermodynamic properties 

In order to facilitate an exact solution, we consider the following infinite-range hopping model which is obtained by 
setting tij in Eq. (1) to a constant properly scaled by the system size tij = t/L so that the thermodynamic limit is 
well-defined. The Hamiltonian is then 



ff = -|E«l«^--z(E«I)(E«^)- (13) 

ij i j 

By substituting the Fourier transform of a^'s defined in Sec. II, bk = J2j e~*"''^aj", one obtains: 

H = -tblbo. (14) 

This model has a very simple spectrum with a ground state with energy —t, and all the other excited states are 
degenerate with zero energy. This particularly simplified spectrum makes an exact solution possible. 

To study the finite temperature properties of this model, we will work with the grand canonical ensemble (GCE), 
in which the chemical potential fi is introduced to control the average density of the system. This model exhibits 
BEG at finite temperature Tq. To determine Tc, we start by considering a system of finite size L, and its occupation 
numbers are: 

{N,^o) = (No) = (Nk) = ^,(_,!) _ ^ for k ^ 0. (15) 

Here (Nq) and {Nk) denote average occupation numbers for the corresponding states in fc-space; P ~ ^ is the inverse 
temperature. From this point on, when we write (-/Vfe), it immediately indicates k 0. The average total particle 
number of the system will be denoted as (N) . To identify Tc, we know in the thermodynamic limit, when T , 
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/i Ek=o = —t, and -^^^ ^ 0. Therefore, (Nk) = ^p}t „ i — -^x^ = (j^) where (n) is the average particle density. So 
we obtain: 

= ln(l + l/(n))- 

Above Tc, { n) — ^ ^ ^-t^ti^i + 7j g,3(-t-^) _2 J from which in the large L limit we can derive that 

M^-Tln(l + ^). (17) 

H has a finite size correction which is negligible above Tc, but will become important below Tc- 
We also know that the partition function of the system in GCE. bears the following form: 

from which it is easy to show that the entropy in GCE. takes the following form: 

S = — ^ — = InZ — — In Z 

dT T dp (19) 

= (1 + (iVo)) ln(l + {Na)) - (No) InNa + {L - 1) [(1 + (iV^) ln(l + (Nk)) - (Nk) In (N^)] . 

Anticipating later relevance, we are particularly interested in the behavior of finite size systems near Tc- For a 
finite system, the chemical potential yU is no longer strictly equal the ground state energy below Tc, but picks up a 
finite size correction Sfi determined by the following condition: 

,-ff/,_l = (20) 

from which we can easily derive that 

(5/x = -Tin ( 1 + — ^ ) . (21) 



(No) 

Consider T = Tc ~ in{i+i/{n)) ' make use of the following fact 



we obtain 



{No) = {N)-Y: (N,) = (A^) - .,ts,Lv (22) 



= - eMt^J- 1 - - (l + L/(A))a'+l/(Ao))-l - ^''^ 



This equation can be solved to give (Aq) as a function of system size T at a given density (n) = (A) /T, at T = T( 



c- 



^^"^ ^ + ^ + ^ ^V{ny + {n). (24) 

Even though this divergent Aq does not affect the thermodynamic behavior of the system, as we will see later it makes 
a (leading) divergent contribution to the mutual information at T = Tc depending on how the system is partitioned, 
or specifically how large is the subsystem size La compared with this VX divergence. 



3.2. Formalism and issues 



In the following part, we will use Peschel's result [32] on the reduced density matrix of a Gaussian state: 



(25) 
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where Gij = {a\aj) is the two point correlation function matrix truncated within the subsystem, and K is the 
normahzation factor. The entropy is given as 

i?A = ^ [( 1 + .9/ ) ln( 1 + .9/ ) - .9i In 9i] , (26) 

where gis are the complete set of eigenvalues of G's (after truncation). Actually this formula also applies to the 
original system. 

We must note that, this formula does not lead to the correct zero temperature limit for the entropy. At zero 
temperature, Gij = {O'l'^j) = ("■)■ Its eigenvalues are all zero except for one: go = {n)L = {N), which gives us a 
non-zero entropy 5t=o = {{N) + 1) \n{{N) + 1) - {N) In {N) = In {N) + {N) ln(l + j^) ~ In {N) at T = 0. This 

reflects the fact that wc arc working with GCE where the particlc-numbcr fluctuation is still permissible at T = 
and the fluctuation amplitude 5N ~ {N). However, as we show below, the mutual information still converges to the 
correct zero temperature limit, the entanglement entropy, at least to the leading order. 
The von Neumann entropy for a subsystem A is given by 

^(GCE) ^ ^ ^^^^^^ + 1)-Na \^Na = I^Na + Na ln(l + -^), (27) 

Na 

where Na ~ {n)LA is the average total particle number in the subsystem A. In the large N limit, the second term 
converges to 1. So the mutual information is given by: 

Em = ^{Ea + Eb~ Sgce r=o) = 3 ( + ^ ) ' (28) 



For Na < Nb, we have 

Eat ^ 

2 

This agrees with Eq. (11) at the leading order. 



Em ^\lnNA + 0{l). (29) 



3.3. Mutual information 

According to our Eq. (25) and Eq. (26), to obtain the von Neumann entropy of the reduced density matrix, 
all what we have to do is to diagonalize the truncated two-point correlation function matrix. Fortunately, within 
this infinite-range hopping model, this is rather simple. For a finite system, we can obtain a general result for all 
temperatures: 

In the above calculation, we have made use of the fact that {Nk) is k- independent. This matrix is easily diagonalized. 
For a system of size L, and a G truncated to a size of La x La denoted by Ga, the eigenvalues are 

51 = + ^^(Nk), 91 = (N,) for / = 2, . . . , La- (31) 

Now the von Neumann entropy of subsystem A can be calculated directly from above result: 

La 

Ea = Y^ ((1 + gi) ln(l + gi) - gi Ingi) 



1=1 



— ^— (32) 



L L ' ' J \ L L 

+ {La - 1) [(1 + (iVfc)) ln(l + (Nk)) - {N^) H{Nk))] 
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Combining the above with Eq. (19), we can obtain the mutual information for a general bipartite system: 



Em ^^{Ea + Eb-S) 



l + ^^ + ^W)ln(l + 



La{No) 



l + ^^^ + ^(A^.))ln(l + 



f LAjNo) Le 
\ L L 

Lb{Nq) , La 



L 



L 



L 



{Nk 



(33) 



Lb{No) , La 



{Nk 



L L 

- (1 + {Nk)) ln(l + {Nk)) + {Nk) In {Nk) - (1 + (^o)) ln(l + (A^o)) + {No) In (TVo) 

Next, we shall discuss the asymptotic behavior of Em in different temperature regions and with different partitions, 
(f) La ^L,T> Tc- 

in this case, {Nq) and {Nk) ~ (n) are both of order one , so ^lA^ ^ o, ii^^ ^ 0, is ^ L. Thus 



Em^O[^ 



(2) La = Lb = k,T > Tc: in this case, Em is reduced to 



(34) 



- {{No) + (7V,))ln((iVo)/2+ (7Vfe)/2) 

-(f + (iV,.))ln(l + (iVfc)) + (7Vfe)ln(7Vfe) ^ ^ 

- (1 + (iVo)) ln(l + (A^o)) + {No) In (A^o) 

which can also be written as an explicit function of L, {n) and T using previous results. 

(3) La ^ L, T < Tc'. in this case, {Nq) becomes a macroscopic number, while {Nk) remains to be of order 1. 
Therefore we isolate the contribution from (A'o), and other terms are of order 0(1): 



EM = lln(l + ^)+Oil). 



If we define (iV^o) = ^^^^ to be the average particle number in the condensate of the subsystem, then 

EM^^ln {NAo)+Oil). 

(4) La — Lb — T < Tc'- the leading contribution is again obtained by keeping {No)^s contribution only: 



Em = - In -(A^o) + 1 + Oil) = - In (A^^o) + 0{l). 



(36) 



(37) 



(38) 



(5) La <^ L, T — Tc'- as we calculated before, (A'^o) diverges as a/L- When La <C L, according to Eqs. (36) and 
(24), we have 



Em = 2 



^ ^ LA^{{n) + l){n) 



Oil) 



■In 



Oil). 



(39) 



For such partition, the scaling behavior of mutual information depends on the ratio If we consider La 

is a small but still finite fraction of L, the scaling behavior of the mutual information still persists: Em = 
iln((n)L^)+0(l). 
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FIG. 1: (Color online) Numerical calculation of mutual information for the infinite-range hopping model for equally partitioned 
systems with average density (n) — 1. The scatters are numerical data, while the dash lines are obtained from our analytic 
results corresponding to the particular temperature. We see exactly what our analytic results tell us: above Tc = the 
mutual information saturates; below Tc, Em ~ I^a; at Tc, Em ~ jLa- Note that the analytic results (dash lines) deviate 
from the numerical calculation because we only keep terms to the subleading order; terms that goes to zero [i.e., of order 
0(-j-!j)] in the thermodynamic limit are neglected. 



(6) La = L/2, T = Tc: by referring to Eqs. (38) and (24), we have 

EM = lln(^\^i{n) + l){n)L^ + 0(1) = 1 ln((n)i^) + 0(1). (40) 

To sum up, from this study, we find that the extensive part of the thermal entropy of the whole system is canceled 
out in the mutual information. Below Tq the mutual information is dominated by contribution from the condensate. 
Even above Tq, contribution from the condensate is of the same order as that from the excited states. The mutual 
information really characterizes the quantum feature of the system. 

To visualize the behavior of mutual information, we present a numerical calculation for the mutual information 
of this model in Fig. (1). This is done by numerically diagonalizing the truncated two-point correlation function 
matrix then computing the von Neumann entropy of the reduced density matrix from those eigenvalues. The system 
is equally partitioned, (n) = 1, so Tc ~ We see exactly what our analytic results tell us: above Tc the mutual 
information saturates; below Tc, Em — \La', at Tc, Em ~ \La- Note that in the plot, the analytic results (dash 
lines) deviate from the numerical calculation because we only keep terms to the subleading order; terms that goes to 
zero [i.e. of order C(-/7j)] in the thermodynamic limit are neglected. 



4. NUMERICAL STUDY OF MUTUAL INFORMATION IN ONE-DIMENSION WITH LONG-RANGE 

HOPPING 



In this section, we shall present our results of numerical study of the mutual information of free bosons living on 
a one dimensional lattice. In this numerical study, we adopt our previous method, and calculate the von Neumann 
entropy of the reduced density matrix from the eigenvalues of the truncated two-point correlation function matrix. 
Throughout this calculation, we hold the average density fixed at (n) = 1 (which means we keep adjusting the chemical 
potential at different temperatures) and consider equal partition only. 

It is well known that for nearest-neighbor (NN) (or other short-range) hopping models whose dispersion relation 
at long- wave length takes the form e(fc) ~ fc^, a finite temperature BEC can only exist in three dimensions (3D). 
However, 3D is in general very challenging for a numerical study that requires large system sizes. Moreover, in 3D the 
mutual information is dominated by area law [10], which renders the logarithmic divergence suggested by our study in 
Sec. Ill sub-leading and thus difficult to isolate. For both of these reasons, it is desirable to study a model in ID with 
BEC at finite T . In ID, the short-range hopping model does not support BEC at finite temperature. To stabilize a 
condensate in ID, we introduce power-law long-range hopping in our free boson model to modify its long-wave length 
dispersion. This is similar to what was done in Ref. [33], in which the authors introduced long-range interaction 
between spins to stabilize magnetic order in ID. The Hamiltonian with long-range hopping is obtained by setting tij 
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in Eq. (1) to the following form, tuned by a parameter 7: 



ij ' -J' k n=l 



L-1 



cos(nfc) 



(41) 



We will show that the long wave-length dispersion is modified to be £^{k) ^ k'' ^ for 7 < 3, as a result if which a 
finite temperature BEC exists for 7 < 2. 

in the thermodynamic limit: 



Consider the eigenenergy function e-y(fc) = — 2t^^^^ cos{nk) 



e^(fc) - -2t J2 



cos(nfc) 



= ~2t Re 



„ink 

Y — 



-2t Re [i^(7, i/c)] , 



where F{'^,v) is the Bose-Einstein integral function [34] defined as: 

1 r 0,7-1 00 

The analytic properties of F{'j,v) near u = are known [34]: 



(42) 



(43) 



r(i-7)i>T-i + 



E 



C(7- 



(7-1)! 



■7-1 , 

E--M«) 



00 

E 

71,^7 — 1 



C(7 " n) 



7^Z, 

7ez, 



(44) 



where C(a;) is the Riemann zeta function. Thus we find that e^(fc) — > k^~-^ for small fc when 1 < 7 < 3. When 7 > 3, 
the low energy dispersion is dominated by the fc^ term. When 7 < 1, s-yik) is not well-defined in the thermodynamic 
limit; in order to have a well-defined thermodynamic limit, the hopping energy t must be properly scaled by the 
system size in this case. 

Next we will consider the thermodynamics of this model with different 7 and demonstrate that for 7 < 2, we 
indeed have a finite temperature BEC. At low temperature, only the small k part of the spectrum is important. For 
1 < 7 < 3 we consider free bosons with a dispersion ak'^~^. Here a = — 2tr(l — 7) is given in Eq. (44). The average 
density of such system in the thermodynamic limit is given by: 



(AO 
L 



1 

2^ 



dk{Nk) = 



1 (7 l-T 

2^ (7 - 1) 



de 



(/3a) — 



- rf- 

2-ie/3£ - 1 27r(7-l) ^7 



-)g 1 (z), 
1 t-i 



(45) 



e''^, /? = i is the inverse temperature, and gv{z) — z-'^e^-i Bose-Einstein integral 



1, which indicates ^-r > 1, 

' 7— 1 ' 



where z — "^'^ « _ 1 

fimction. To have a finite temperature BEC, {N)/L = (n) must remain finite when z 
because for v < 1, 5^(1) diverges. 

To have a better understanding of the thermodynamics of this model, in Fig. (2) we present a numerical calculation 
of Tc- The exact average density of the system is: 



(AO 
L 



1 

2^ 



dk- 



,/3(e^(fc)-M) - I' 



(46) 



where e~f{k) is the eigenenergy function given in Eq. (42). Tq is computed by setting /i = e-y(O) and then solving this 
equation numerically. According to Fig. 2, Tq grows monotonically from to cx3 as 7 goes from 2 to 1. The divergent 
behavior of Tc as 7 ^ 1 is a consequence of the divergent bandwidth in that limit. 

According to our study of the infinite range hopping model, above Tc the mutual information should saturate as the 
system size grows. Below Tc, the mutual information has a scaling behavior Em ~ ^ InL^; for T = Tc Em — \ In^A 
for equally partitioned system in that model. We expect the In La scaling behavior both below Tc and at Tc to persist 
in the long-range hopping model. As we shall see later, this is indeed the case. However, the details of the scaling 
behavior (i.e. the prefactor) can be different for different 7. To study this scaling behavior, we fix the temperature and 
examine the mutual information as a function of system size. This is desirable because, if our conjecture according 
to the study of infinite-range hopping model is correct, the mutual information will be proportional to Ini^ when 
T<Tc. 
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FIG. 2: (Color online) Numerical calculation of Tc (line + symbols) for the long-range hopping model in the thermodynamic 
limit. Tc is measured in unit of the hopping energy t which is set to 1. As one can see, Tc grows monotonically from to cxd 
as 7 goes from 2 to 1. The divergent behavior of Tc as 7 —> 1 is a consequence of the divergent bandwidth. 
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FIG. 3: (Color online) Mutual information of the nearest neighbor hopping model plotted against subsystem size on a logarithmic 
scale. Average density is set to (n) = 1, and the system is equally partitioned, L — 2La- The black dash line is Em ~ § IuLa. 
This line will be in other graphs for comparison as well. Clearly, the mutual information saturates when the system size grows 
large enough. 



Before we demonstrate our results in the long-range hopping model, first let us verify our analysis in the NN 
hopping model in which no BEC would occur. This case actually corresponds to 7 ^ 00. The Hamiltonian is given 
by restricting the hopping in Eq. (1) to the nearest neighbors only: 

H^.n = -tY,4cj- (47) 

iv) 

Figure 3 is a linear-log plot of mutual information against subsystem size at different temperatures. Throughout our 
study we shall consider equal partition only, i.e.. La = L/2. The average density is also set to (n) = 1 here, so are 
the other results we will show later. Clearly, at a fixed temperature, the mutual information saturates as the system 
size grows. At low temperatures, small systems can be considered in the zero temperature limit. This leads to the 
mutual information growing as ~ ^ IuLa, until saturation kicks in. 

Next we consider 1 < 7 < 2. Now the finite temperature transition emerges, and the signature for the transition 
in mutual information - the logarithmic scaling with (sub-)system size also emerges. In Fig. 4, we plot the mutual 
information (scatters) for 7 = 1.7 at different temperatures. The best fit (cyan dash line corresponding to the diamond 
data points) for /3 = 0.5 > f3c gives Em = 0.2405 Ini^ -I- 0.214. At (Jc = 0.297, the scaling behavior is fit (magenta 
dash line) as Em = 0.1226 InL^ -I- 0.1688. Both behaviors agrees qualitatively with what has been suggested by our 
analytic study of the infinite-range hopping model. When the temperature is well below Tc, we have the logarithmic 
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FIG. 4: (Color online) Mutual information of the long-range hopping model with the parameter 7 = 1.7 as a function of 
subsystem size on a logarithmic scale, at various (inverse) temperatures. The average boson density (n) is set to 1, and the 
system is equally partitioned, La = L/2. The scaling behavior for inverse temperature /3 = 0.5 goes as Em — 0.2405 In La+0.214 
(cyan dash line corresponding to the diamond data points). At the transition point, /3 — Pc = 0.297, the scaling law is fit to 
be Em = 0.1226 In La + 0.1688 (red dash line corresponding to the right triangle data points). 



scaling behavior: a set of parallel linear lines on this logarithm scale plot for different temperatures. However, the 
prefactor is significantly different from that of the infinite range hopping model which is ^. In fact, by calculating 
mutual information of different 7's, we find that the prefactor varies as 7 changes. For very low temperature, small 
systems again effectively fall into the zero temperature region, and the mutual information restores to the ~ ^ In La 
behavior as in zero temperature. But when the system size becomes large it crosses back to the finite temperature 
scaling behavior again. This is evident for /3 > 1 in the figure. For temperature close to but still below Tp, small 
systems behave differently: the mutual information scales more like the line for /3 = j3c-, and it bends up as the system 
size increases and finally crosses back to its genuine behavior below Tc- For temperature close to Tc but now above, 
small systems behave the other way: the mutual information bends downwards and saturates at large system size. 
j3c serves as a very distinctive boundary between the above two different bending behaviors. The latter two bending 
features are also present in Fig. 1, our numerical verification of the infinite-range hopping model. But the first feature 
at very low temperature is missing in Fig. 1 since in that case the mutual information scales the same way as the 
entanglement entropy. 

Very similar behaviors were observed for the entire range 1 < 7 < 2. Representative results are presented in Figs. 
5 and 6 for 7 = 1.5 and 1.3 respectively. We thus conclude that for the entire range 1 < 7 < 2, mutual information 
saturates for T > Tc, while it diverges logarithmically with increasing subsystems size, for both T < Tc and T = Tc- 
The coefficients in front of the logarithms are 7-dependent. 



5. SUMMARY AND CONCLUDING REMARKS 



In this paper we have studied entanglement properties of free non-relativistic Bose gases. At zero temperature, all 
particles fall into the ground state, and we find the entanglement entropy diverges as the logarithm of the particle 
number in the subsystem. At finite temperatures, we studied the natural generalization of entanglement entropy - 
the mutual information. We find the mutual information has a similar divergence in the presence of a Bose-Einstein 
condensate. When the system is above Tc or does not have a condensate, the mutual information saturates for 
large subsystem size. It should be noted that for the special models we studied in this paper there is no area-law 
contribution to the mutual information, thus the contribution from the condensate, when present, dominates the 
mutual information. In more generic models in two- or three-dimensions where an area-law contribution is present, 
we expect such logarithmic divergent contribution from the condensate to be present as a sub-leading term in the 
subsystem-size dependence of the entanglement entropy and mutual information. 

Physically it is easy to understand why the condensate makes such an important contribution to entanglement. 
First of all, BEC is intrinsically a quantum process, just like entanglement reflects the intrinsically quantum nature of 
the system. More specifically, when a (macroscopically) large number of particles occupy the same state (at k = 0), 
they are necessarily delocalizcd throughout the sample, giving rise to entanglement between blocks. 

Just like in our previous work on a very different system [27], our results here suggest that conventional ordering. 
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FIG. 5: (Color online) Mutual information of the long-range hopping model with the parameter 7 = 1.5 as a function of 
subsystem size on a logarithmic scale, at various (inverse) temperatures. The average boson density (n) is set to 1, and the 
system is equally partitioned, La ~ L/2. The mutual information for /3 = 1 is fit to scale as Em — 0.324 In L -I- 0.445 (orange 
dash line). At /3c ~ 0.16843, we observer a weaker scaling behavior which is fit to be Em — 0.064 In La + 0.084 (red dash line). 
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FIG. 6: (Color online) Mutual information of the long-range hopping model with the parameter 7 = 1.3 as a function of 
subsystem size on a logarithmic scale, at various (inverse) temperatures. The average boson density (n) is set to 1, and the 
system is equally partitioned. La = L/2. The best fitting for /? — 0.5 > /3c (cyan dash line) gives Em = 0.378 In La + 0.2. At 
(3c = 0.0954, the scaling behavior is fit (red dash line) as Em = 0.023 In La + 0.03535. 



like BEC, makes a logarithmic contribution to entanglement. One thus needs to take caution when using entanglement 
as a diagnostic for exotic phases (such as topological phases) or quantum criticality. 
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